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Ab initio calculations using the Perturbed Ion model, with correlation contributions included, are 
presented for nonstoichiometric (NaI)„Na"'" and (CsI)„Cs"'' (n<14) cluster ions. The ground state 
and several low-lying isomers are identified and described. Rocksalt ground states are common and 
appear at cluster sizes lower than in the corresponding neutral systems. The most salient features of 
the measured mobilities seem to be explained by arguments related to the changes of the compactness 
of the clusters as a function of size. The stability of the cluster ions against evaporation of a single 
alkali halide molecule shows variations that explain the enhanced stabilities found experimentally for 
cluster sizes n=4, 6, 9, and 13. Finally, the ionization energies and the orbital eigenvalue spectrum 
of two (NaI)i3Na^ isomers are calculated and shown to be a fingerprint of the structure. 

PACS numbers: 31.20.Gm 36.40.+d 61.46.+W 



I. INTRODUCTION 

Alkali halide clusters have received substantial experi- 
mental and theoretical attention. From the experimental 
side, they are relatively easy to form, and their ionic 
bonding character has allowed the development of sim- 
ple models for the interaction between the ions that form 
the cluster. They provide a nice opportunity to study 
the emergence of condensed matter properties. An im- 
portant point in this respect is the determination of the 
most stable isomers for each cluster size. 

In this paper, we present the results of ab initio cal- 
culations of the structures and stabilities of (NaI)„Na''" 
and (CsI)„Cs+ clusters with n=l-14. Large clusters of 
these materials were first produced by bombardment of 
crystalline targets with Xe+ ions ancLdetected by means 
of secondary-ion mass spectrometry.EJu The mass spec- 
tra showed anomalies in the cluster intensities for certain 
values of n, which were tentatively interpreted as reveal- 
ing the enhanced stability of "cuboid-like" structures. 
The results were explained in terms of a direct emis- 
sion model, which assumes that "cuboid-like" cluster ions 
are directly sputtered from the crystal. Ens et alu per- 
formed time-of-flight measurements for (CsI)„Cs+ clus- 
ters, again produced by bombardment of Csl crystals, 
and considered different observation times after cluster 
emission. In that way, they found that the anomalies in 
the mass spectra were a consequence of evaporative de- 
cay a/ier production of the cluster ions, and that the pre- 
ferred decay channel was the ejection of a Csl molecule. 
Anomalies in the mass spectra of clusters formed by 
the inert-gas -condensation technique were observed by 
Pflaum et al.a. Twu et al. published mass spectra of 
sodium chloride, sodiurp, iodide, cesium chloride, and ce- 
sium iodide cluster ionsci produced by laser vaporization. 
They observed the same magic numbers for all those ma- 



terials, namely n=4, 6, 9, and 13 in the small size range. 
The differences in the detailed structure of the mass spec- 
tra were attributed to the relative sizes of the ions making 
up the clusters. All those techniques, however, give no 
direct information of the cluster shapes. Drift cell experi- 
ments, which measure the mobility of the cluster ions and 
stand as a promising technique for structural analysis of 
clusters, have been performed to studyr-the structures of 
small covalent and metallic clusters. Ert2l More recently, 
ion mobility CKoeriments have been performed for alkali 
halide .chisteralj and have revealed isomerization transi- 
tionslay in (NaCl)„Cl-. j-. 

It was pointed out by MartinlU that precise measure- 
ments of the photoionization spectrum should also help 
in determining the structure of clusters of ionic materi- 
als. The reason is that, due to the strong localization 
of the electrons in closed-shell alkali halide clusters, the 
structure of the photoionization spectrum should be a 
fingerprint of the structure of the cluster, giving infor- 
mation on the set of inequivalent anions. Li et al. used 
this idea to correlate the optical absorption spectra of 
(CsI)„Cs+ with ±hc_structures obtained by using a pair 
potential model.EjtZl Photoelectron and photoionization 
spectroscopy as well as theoretical studies have also been 
performed to study the relation between the cluster struc- 
ture and the localization ipDde.for small alkali halide clus- 
ters with excess electrons ,E3"EiJ and to study the structure 
and emergence of metallic propertieaJji alkali-rich alkali 
halide and alkali hydride clusters. EHtEj 

Some theoretical calculations for alkali halide clus- 
ters have bpcip-.based on phenomenological pair poten- 
tial modelsJlJ'EZlEZrc^ Such simplified models have been 
successful in explaining their main characteristics, and 
furthermore, are very useful to find the different local 
minima in the potential energy surfajce._j4& initio calcula- 
tions performed by Ochsenfeld et aZ.ErE3 used molecular 
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orbital methods including correlation at the MP2 levd_ 
We have used the ab initio perturbed ion (PI) model,Ej 
which is formulated within the restricted Hartree-Fock 
(RHF) approxinaation, in studies of neutral stoichiomet- 
ric alkali halidcEj £3 and (MgO)„ clusters, lZI and in a 
prcliiainary study of nonstoichiometric (NaCl)„Na+ clus- 
ters.E3 In some of those studies correlation contributions 
were included using an unrelaxed CouLsmb-Hartree-Fock 
(uCHF) model proposed by Clementi.OEJ The PI model 
represents a major advance with respect to pair poten- 
tial methods and provides an alternative description to 
the molecular orbital methods. In this paper we present 
an extensive study of charged (NaI)„Na~'" and (CsI)„Cs~'' 
clusters with n up to 14. The results are aimed to assist 
in the interpretation of the experimental investigations of 
the structure of alkali halide clusters, as provided by the 
ion mobility studies, or of the relation between structure 
and photoionization spectrum. 

The remainder of this paper is organized as follows: 
In Sec. ^ we briefly review the computational method 
used in this study. In Sec. Ill we report the results for 
the isomer structures, relative stabilities, and ionization 
potentials. Section |^ summarizes our conclusions. 



II. COMPUTATIONAL METHOD 

The ab initio perturbed ion model,^! as adapted to 
the study of ionic clusiep. has been described at length 
in our previous works EjEj In brief, the HF equations of 
the cluster are selfconsistently solved in localized Fock 
spaces, by breaking the cluster wave function into local 
nearly orthogonal ionic wave functions and accounting 
for ion-cluster consistency. The average binding energy 
per ion of the (AX)„A"'" cluster is given by 
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[nEoiX-) + {n + l)EoiA+) - £;ciustcr], 

(1) 



where Eo(A+) and Eo(X~) are the energies of the free 
cation A+ and anion X~ , respectively. The localized na- 
ture of the PI model wave functions has some advantages. 
In weakly overlapping systems, the correlation energy is 
almost intra-atomic. In this paper, the correlation en- 
ergy correction is obtained through dementi's Coulomb- 
Hartree-Fock method.E^cJ The PI model also allows for 
the development of efficient computational codesuiLwhich 
make use of large multi-zeta Slater-type basis set^o for 
the description of the ions. Our calculations have been 
performed using the following basis sets: (Ils9p5d) for 
Cs+, taken from Ref. |3|; (5s4p) for Na+, and (Ils9p5d) 
for I^, both taken from Ref. As input geometries for 
the optimization of the atomic structure we have consid- 
ered thp .stj.uctures obtained from pair potential calcu- 
lations.oEZlca'CJ Those input geometries have been fully 
relaxed for (NaI)„Na+ (n < 6) and (CsI)„Cs+ (n < 4) 
clusters, that is, the total cluster energy has been min- 
imized with respect to variations in all the (3N-6) inde- 
pendent coordinates, where N is the number of ions. For 
larger clusters, full relaxations are expensive and we have 



relaxed these structures with respect to a limited num- 
ber of relevant parameters, which depend on the spatial 
symmetry of each isomer. A downhill simplex algorithm 
has been used in all the optimizations .C3 All the cluster 
energies are converged within a precision of 1 meV. 



III. RESULTS 

A. Structures of Isomers 

The calculated structures of small (NaI)„Na+ and 
(CsI)„Cs+ clusters are shown in Figs. 1 and 2, respec- 
tively. Small spheres represent cations and large spheres 
represent anions. The most stable structure, or ground 
state (GS) is shown on the left side for each cluster size. 
The other structures represent low-lying isomers. Below 
each isomer, the energy difference with respect to the 
ground state is given. The structures exhibit several dis- 
tinct motives which can be roughly classified in chains, 
rocksalt pieces, and rings (mainly hexagonal). It is pos- 
sible that other isomers could exist between the ground 
state and the low-lying isomers plotted in the figure, 
since our search has been limited mainly tritho -structures 
provided by pair potential calculations,oEZlE3o and the 
possibility of overlooking some isomers can not be ex- 
cluded. For n > 10, the rocksalt motives consistently 
dominate the ground state and the crystalline structure 
emerges, although signs of the appearance of the rock- 
salt structure are also found for some of the clusters with 
n < 10. Let us describe first the structures obtained for 
the smaller cluster sizes. A linear chain is obtained as 
the ground state of (NaI)2Na"'". Bent chains, a rhombus 
with a cation attached, and a three-dimensional isomer 
are found less stable. All those structures are closer in 
energy for (CsI)2Cs''" and its three-dimensional GS is de- 
generate with the linear chain. The same GS structure is 
found for both materials at n=3, namely a cube with an 
anion removed. The linear chain is still the second iso- 
mer for Nal, while it is found at a higher energy in Csl. A 
planar rocksalt piece plus a cation appears as a high lying 
isomer. The rest of the isomers are obtained by attaching 
a molecule in several ways to some of the n=2 structures. 
The n=4 GS can be described as a quasi-two-dimensional 
sheet, which is quite curved in Nal but almost planar for 
Csl. The chain isomers are not competitive anymore. A 
cube with a cation attached is obtained as the second 
isomer. Again an essentially two-dimensional sheet, de- 
rived from the GS structure of (NaI)4Na"'" is the GS of 
(NaI)5Na+. This isomer is still more stable than the 
three-dimensional (3D) structures, although the 2x2x3 
rocksalt piece with a corner anion removed is energeti- 
cally very close. This notation indicates the number of 
ions along the x, y and z directions, respectively. For 
Csl, a cube with a linear chain on-top is obtained as 
the GS. The tendency for three-dimensional structures is 
stronger for the Cs-clusters. Rounded cages with quad- 
rangular and hexagonal faces form the GS for n=6 and 
n=7 in Fig. 1. Rocksalt isomers are still high in energy, 
although an hexagonal isomer is only 0.02 eV above the 
GS in (NaI)6Na+. GS structures with higher coordina- 
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tion are obtained for Csl at these cluster sizes: a centered 
hexagonal prism for n=6, and a complex structure con- 
taining a cube for n=7. We can notice that the centered 
hexagonal prysm is also related to the rock-salt struc- 
ture since the central cation is coordinated to six ions. 
Both in the ground state and in the low-lying isomers 
discussed until now (n < 7) we appreciate a tendency 
to distorted structures that, in the more extreme cases 
leads to elongated clusters. This is driven by the excess 
positive charge. The distortion (elongation) lowers the 
repulsion between cations. 

A 3 X 3 X 2-like structure with an anion missing from the 
center of a basal plane is obtained as the GS for n=8. The 
absence of this cation induces a distortion of that basis to 
an octogonal ring. A structure with an eight-coordinated 
cation is close to the GS for Csl. Eight is the coordina- 
tion number of solid Csl. Structures obtained by adding a 
cation in several ways to the hexagonal prismatic form of 
(Nal)9 (which is the GS of the neutral clusterEj) are the 
most stable (NaI)9Na"'" isomers. In the GS, this cation 
induces a strong cluster deformation. The (CsI)9Cs"'" GS 
is also related to an hexagonal prysm. Cuboid-like rock- 
salt structures are less stable in both systems, but we 
can again notice that the upper half of the hexagonal 
prysm is distorted in (CsI)9Cs~'", and contains an inner 
cation with coordination six (the rock-salt coordination) . 
Starting with n=10 the emergence of cuboid- like rock- 
salt features becomes apparent. From n=10 to n=14, 
all the Nal GS clusters have the rocksalt structure. No- 
tice in particular the high stability for n=13. More open 
structures with lower average coordination, or based on a 
stacking of rings are less stable. The same can be said of 
Csl, with due exception of n=12, where a hexagonal pris- 
matic structure (although with a central six-coordinated 
cation) is more stable. 

The structural trends of neutral alkali halide clusters 
have been studied in Ref. Those trends were ratio- 
nalized in terms of the relative ion sizes. As the ratio 
rc/rA between the cation and anion radii increases, the 
tendency to form rocksalt fragments becomes enhanced. 
Nal has a small value of this ratio (0.44) while Csl has a 
large ratio (0.76). In spite of this difference, the charged 
nonstoichiometric clusters behave in a similar way. Rock- 
salt pieces appear early; about three quarters of the clus- 
ters between n=3 and n=14 have a rocksalt fragment, or 
a closely related structure, as the GS. The exceptions oc- 
cur for n=6, 7, 9; in these cases the rocksalt isomers have 
one or two low-coordinated cations. One can notice the 
influence of nonstoichiometry and net charge: the per- 
cent of GS rocksalt structures in neutral (Nal)„ is only 
about one half.Ea The reason seems to be that hexagonal 
prysmatic structures are less competitive for the charged 
nonstoichiometric clusters. A perfect prysm formed by 
hexagonal rings has equal number of cations and anions, 
so only defect-structures, obtained by removing an anion 
or by adding a cation, can be built for nonstoichiomet- 
ric clusters. We find one example of the first type in 
one of the isomers of (NaI)5Na"'". A cation can be added 
on top of a terminal ring or in the interstitial hole be- 
tween two hexagonal rings: isomers of both types exist 
for (NaI)6Na+. These only become competitive when the 



rocksalt isomers are very unstable, like for n=6, 9. 

Pair potential calculations have beea, performed by 
Diefenbach and Martinet and Li et alt3 for these two 
systems. A comparison with the PI results shows 
some discrepancies. The pair potential calculations for 
(NaI)„Na+ with n=8, 10, 11, 12, and 14 predict rather 
complex GS structures (independently of the use or not 
of polarization terms) which correspond to some of the 
isomers in fig. 1, whereas the PI calculations lead to 
rock-salt-type structures. For (CsI)„Cs"'" the discrepan- 
cies are rather insignificant. In the PI niei|del the bind- 
ing energy of the cluster can be writtenE3 as a sum of 
classical and quantum interaction energies between the 
ions plus a term that accounts for the radial deforma- 
tion of the electronic cloud of the free ions (in practice 
the anions) in response to the environment. These energy 
contributions contain quantum many electron terms that 
in principle describe the interatomic interactions better 
than the empirical potentials. There is also a second 
type of discrepancies between the PI and pair potential 
calculations, although much less significant. These occur 
when the rigid-ion and polarizable-ion model potentials 
disagree with each other and the PI model agrees with 
one of them. These cases are found for (NaI)„Na"'" with 
n=2, 5, and for (CsI)„Cs+ with n=2, 7, 8, 11, 12. In 
those cases the PI model agrees sometimes with the rigid 
ion model and sometimes with the polarizablc ion model, 
but the two isomers involved are generally close in energy, 
so the nonuniformity of the agreement can be ascribed 
to the small energy differences involved. In addition one 
has to bear in mind that, due to one basic assumption 
of the PI model (spherically symmetric electron density 
clouds, centered on the nuclei), the polarization contri- 
bution coming from dipolar terms is absent in this model. 

Kreisle and coworkers^ have studied the mobility of 
(NaI)„Na"'" and (CsI)„Cs^ clusters under the influence 
of an electric fleld in a chamber filled with helium gas. In 
these experiments the mobility is larger the lower the He 
scattering cross section by the cluster, and this cross sec- 
tion is inversely related to the compactness of the cluster. 
Kreisle and coworkers have plotted the inverse mobilities 
as a function of cluster size. Some salient features are 
common to the two curves and, in our opinion, can be 
related to the structural features found in Figs. 1 and 2. 
The main feature is a clear drop in the inverse mobility 
between n=12 and n=13. In fact, the inverse mobility be- 
comes a local minimum for (NaI)i3Na^. It is suggestive 
to associate the high mobility of n=13 to its compact 
"perfect cube" form. Other feature is a visible change 
(decrease) of the slope of the inverse mobility curve at 
n=4. It is tentative to associate this to the change from 
the two-dimensional to more compact three-dimansional 
character of the ground state of (CsI)„Cs"'" between n=4 
and n=5. Although the calculated GS of (NaI)5Na+ is 
planar, there is a low-lying isomer, only 0.02 eV higher 
in energy, that could easily be present in the beam and 
contribute to increase the mobility. 
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B. Relative Stabilities as a Function of Cluster Size 
and Comparison with Experiment 

Tijo-|experimental mass spectra of alkali halide cluster 
iona^lia show intensity anomalies which reflect the special 
stability of some cluster sizes. In order to study the rel- 
ative stability of (NaI)„Na"'" and (CsI)„Cs"'" cluster ions, 
we plot in Fig. 3 the average binding energies per ion (eq. 
1) as a function of n. Maxima or pronounced changes of 
slope in these curves are considered as indicators of en- 
hanced stability. Clear maxima at n=4 and n=13, and 
main slope changes at n=6 and n=10, are obtained for 
(NaI)„Na"''. For (CsI)„Cs~'" a maximum is apparent at 
n=13, and main slope changes occur at n=4, 6, and 9. 
Those features correlate with the observed abundaBCie 
maxima. The most prominent observed maximumLTEl 
is n=13. The magic numbers at n=4, 6 and 9, and 
specifically the enhanced abundances of (NaI)6Na+ and 
(NaI)9Na+ clustersjj are less pronounced. The only dis- 
crepancy between experiment and theoretical predictions 
occurs for (NaI)9Na"'". However, the slope changes in 
Fig. 3 are so weak that the average binding energies are 
not the best indicators of the enhanced stability of some 
magic clusters. 

The quantity Ehind measures the cluster stability with 
respect to the infinitely separated ions. The experiments 
indicate, however, that the abundance mass spectrum 
should be probably best explained in terms of the st 
ity against evaporation of an alkali halide molecule.! 
The energy required to remove a molecule AI from an 
(AI)„A+ cluster ion (A = Na, Cs) is given by 



^evaporation — 

-E,ius[{AI)nA+] 
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The evaporation energies are plotted in Fig. 4. A sharp 
increase in the evaporation energy between n=14 and 
n=13, between n=10 and n=9, between n=7 and n=6, 
and finally between n=5 and n=4 is evident for the 
(CsI)„Cs+ clusters. This means that evaporative cooling 
will result in enrichment of clusters with n=4, 6, 9 and 13 
in the beam. The results are similar for (NaI)„Na"'', pre- 
dicting enrichment of clusters with n=4, 6 and 13, but a 
discrepancy with respect to experiment is again obtained 
since enrichment is predicted for n=10, instead of n=9. 
In an attempt to resolve this discrepancy we have also 
plotted in Fig. 4 a "vertical" evaporation energy. This 
is defined by removing from the parent cluster (size n) 
the least-bound molecule (this can be identified in the PI 
model, since the total binding energy of the cluster can 
be separated into a sum of ion contributions; see Refs. 
|33|-|36| for details), and relaxing the resulting structure 
(size n-1) to its nearest local minimum. This is in many 
cases not the ground state of the (n-l)-cluster and the 
difference between the "adiabatic" and "vertical" evapo- 
ration energies in Fig. 4 accounts for this fact. In spite 
of this difference the use of vertical evaporation energies 
leads to the same predictions for (Cs^nCs"*" as before, but 
changes the predictions for (NaI)„Na"'" to improve agree- 
ment with experiment for n=9, which is a maximum in 
the curve of the vertical evaporation energies. The in- 
terpretation is that, although the adiabatic evaporation 



of a molecule from (NaI)ioNa+ costs more energy than 
adiabatic evaporation from (NaI)9Na"'", there are in both 
cases isomeric forms of the (n-l)-clusters with i) a struc- 
ture similar to that of the n-cluster and ii) large energy 
barriers between those isomeric forms and the ground 
state structure of the (n-l)-cluster, such that the vertical 
evaporation from (NaI)9Na^ is larger. In summary, our 
calculations suggest the possible relevance of isomers of 
the (n-l)-cluster with a structure similar to that of the 
GS of the n-cluster to explain the details of the mass 
spectra when evaporative cooling is involved. This point 
deserves further investigation. The main magic numbers 
n=4 and n=13 are a consequence of the enhanced stabil- 
ity of very symmetrical rocksalt structures, in two and 
three dimensions, respectively. On the other hand, n=6 
and n=9 are "fine structure" peaks of the spectra and the 
explanation in terms of structural features is less evident. 
These occur for (CsI)„Cs+ because structures are formed 
that optimize the value of the Madelung energy more effi- 
ciently than for neighboring cluster sizes. (CsI)9Cs"'" has 
some highly coordinated ions: one cation with coordina- 
tion 6 and three anions with coordination 5. (CsI)6Cs^ 
also contains one cation with coordination 6. At the 
same time the lowest coordination found in these two 
clusters is 3. In contrast, some of the neighbor clusters, 
like (CsI)5Cs+ and (CsI)ioCs"'" contain some ions with 
coordination 2. Those highly coordinated structures are 
less competitive for (NaI)6Na+ and (Nal)gNa"''. Figure 
3(a) also shows, for comparison, the binding jEnergies per 
ion of neutral stoichiometric (Nal)„ clusters.L3 The local 
maxima occur for n=6, 9 and 12, and have been associ- 
ated to the formation of compact structures with large 
atomic coordination compared to clusters with n-l-1 and 
n-1 molecules. 



C. Ionization Energies and Structure 

In previous studies of alkali-halide clusters ^'ll we 
have analyzed the variation of the ionization potential 
(IP) with the cluster size. The vertical IP was calculated 
in the Koopmans' approximation as the binding energy 
of the lowest bound electron in the cluster, which is of 
course located on a specific anion. Here we investigate 
the relation between the geometrical structure and the 
spectrum of electronic states for different isomers of the 
same cluster. This could provide a way to distinguish, 
isomers, already explored for other types of clusters.^ 
In the PI model the electrons near the HOMO level are 
localized on distinct anions, and the different eigenval- 
ues arise from the different atomic environment around 
nonequivalent anions. Thus, the set of one-electron en- 
ergy eigenvalues characteristic of each isomer cap_be con- 
sidered as a fingerprint of its structural shape.Eil As an 
example, we present in Fig. 5 the orbital energy spectrum 
corresponding to the two isomers of (NaI)i3Na"'" given in 
Fig. 1. It is apparent that the two isomers have quite 
different spectra. These could be measured by photoelec- 
tron spectroscopy, and in principle it could be possible to 
determine the structure of the different isomers present 
in a mass-selected beam by comparing the experimental 
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spectra with theoretical results. In our example, only 
two peaks are expected in the lowest energy part of the 
spectrum for the case of the GS isomer, because the sym- 
metry of this structure induces high degeneracies. One of 
the peaks corresponds to removing one electron from any 
of the twelve surface anions (labelled A in the figure) and 
is identified with the vertical IP. The other corresponds 
to the removal of one electron from the central anion (la- 
belled B) with coordination six. The second isomer has 
a larger number of inequivalent anions and the spectrum 
is broader. Besides, the ionization energy is 1.3 eV lower 
than for the GS isomer. Similar "fingerprints" distin- 
guish different isomers for other cluster sizes. 

IV. SUMMARY 

The determination of the structures of alkali halide 
cluster isomers is a challenging, subject for present- 
day experimental techniques .EilO Theoretical calcula- 
tions can throw light on these problems. In this pa- 
per, we present ab initio calculations of the structures 
and stabilities of (NaI)„Na+ and (CsI)nCs+ cluster ions 
with (n<14). Starting from several local minima found 
with phenomenological pair potential models, we have 
used the Perturbed Ion Model (with correlation included) 
in order to determine the ground state and some low- 
lying isomers. Our results indicate an early formation 
of rocksalt fragments. The rocksalt features appear 
at values, [Of n lower than in the corresponding neutral 
species.tJ'Ea Arguments related to the compactness of 
some clusters appear to be able to explain the main fea- 
tures obtained lisi the mobility measurements of Kreisle 
and coworkers.y The nw-^ spectra obtained by several 
experimental techniquesLTO show evidence of enhanced 
population for cluster sizes n= 4, 6, 9, and 13. Our cal- 
culations confirm the enhanced stability of these clusters, 
namely, those clusters are very stable against evapora- 
tion of a single molecule. We have investigated the pos- 
sibility to determine isomeric structures by comparing 
experimental photoelectron spectra with those obtained 
theoretically. As an example we have shown how the 
spectra of orbital energy eigenvalues of two (NaI)i3Na+ 
isomers depend on the structure. With all these results 
in mind, we are confident that the identification of the 
isomer structures of clusters of ionic materials can be 
feasible in a near future if experimental and theoretical 
efforts work together. 



Captions of figures 

Figure 1. Lowest-energy structure and low- lying iso- 
mers of (NaI)„Na+. The energy difference (in eV) with 
respect to the most stable structure is given below the 
corresponding isomers. 

Figure 2. Lowest-energy structure and low-lying iso- 
mers of (CsI)„Cs"'". The energy difference (in eV) with 
respect to the most stable structure is given below the 
corresponding isomers. 

Figure 3. Binding energy per ion as a function of the 
cluster size of (NaI)„Na+ (circles) and (Nal)„ (squares) 
in panel (a) and of (CsI)„Cs+ in (b). 

Figure 4. Adiabatic (circles) and vertical (squares) 
energies required to evaporate a neutral molecule from 
(NaI)„Na+ (a) and (CsI)„Cs+ (b) clusters as a function 
of n. 

Figure 5. Spectrum of orbital energy eigenvalues for 
the two isomers of (NaI)i3Na+ studied. Inequivalent an- 
ions are labelled with different letters. The vertical scale 
gives the number of equivalent anions of each type. 
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